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Abstract 

We identify models of supersymmetric hybrid inflation in which the tensor- 
to-scalar ratio, a canonical measure of gravity waves produced during infla- 
tion, can be as large as 0.03 or so, which will be tested by the Planck satel- 
lite experiment. The scalar spectral index lies within the WMAP one sigma 
bounds, while \dn s / d\nk\ < 0.01. 



Introduction 

Supersymmetric (SUSY) models of hybrid inflation jl|, 0, IE 0, 0, ls| are 
enticing due to their ties to grand unified theories (GUTs) and mainstream 
particle physics. In addition to utilizing the pervasive framework of SUSY, 
these models naturally incorporate the breaking of a gauge group G into 
the inflationary mechanism, whether at the end of inflation (standard sce- 
nario) or during inflation (shifted scenario 0, Q). The gauge symmetry G 
may be associated with a GUT, such as flipped SU(5) [9|, or perhaps an 
extension of the Standard Model (e.g. an extra U(1)#-l symmetry). With 
the recent launch of the Planck satellite, and in anticipation of the new re- 
sults it promises to yield, a thorough exploration of the predictions of such 
highly motivated models may bring us closer to understanding the nature of 
inflation. 

Primordial gravitational waves in the inflationary epoch can source ten- 
sor fluctuations in the cosmic microwave background (CMB). A convenient 
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parametrization of these tensor modes is obtained by comparison to the scalar 
perturbations, 

A 2 

(1) 



A* 

This 'tensor-to-scalar ratio' is observably large only if the inflaton field 
changes amplitude over an interval near to or greater than the Planck scale, 
owing to the Lyth bound 



r<0.0 8 (^) , (2) 
\m P J 

where nip ?a 2.4 x 10 18 GeV is the reduced Planck mass. This constraint most 
readily allows sizable r- values if A0 3> mp. Recently, some progress has been 
made in devising realistic inflationary models having observable r (~ 0.01 



or so) in which the inflaton amplitude is sub-Planckian [12j. Indeed, we will 
require that the field amplitude remain below the Planck scale in order to 
ensure that supergravity corrections remain under control. 

Thus far, models of SUSY hybrid inflation have predicted relatively tiny 
values of r (< 10 -4 ), leaving little hope of observation by current or future 
experimental endeavors (see, for example, Refs. QHHEDi)- In this letter, 
we will show that the inclusion of higher order corrections from multiple 
sources can open up a region of parameter space that can support r-values 
large enough to be observed by the current Planck satellite observatory. 

Supersymmetric Hybrid Inflation 

Within the framework of supersymmetry, hybrid inflation is achieved via 
the superpotential [H 0] 

W = kS($$ - M 2 ) , (3) 

where S is a gauge singlet superfield whose scalar component acts as the 
inflaton, and $, $ are gauge conjugate superfields transforming nontrivially 
under some gauge group G. It is often desirable to include an additional U(l) 
'.R-symmetry' in SUSY theories (e.g. to curb runaway proton decay). In this 
case, S transforms under U(1)k in the same way as W, and Eq. (J3]) is the most 
general superpotential consistent with U(l)# and G at the renormalizable 
level. 
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We are interested in exploring the possibility of large primordial gravity 
wave amplitudes produced in these models of inflation. Consequently, it 
is insufficient to consider only global SUSY; supergravity (SUGRA) effects 
must be taken into account. The Kahler potential may be written as an 
expansion in powers of 1/mp: 

lrfl2 1*12 1*12 K S |S1 4 K$ |$| 4 ^l^l 4 

K= \S 2 + $ 2 + $ +-r L r + T L r + T L T 

4 trip 4 mp 4 mp 

\S\ 2 \<&\ 2 \S\ 2 \$\ 2 |$| 2 |$| 2 k ss \S\« 

+ Ks$ 5 h 5 h 5 1 t — h • • ■ , (4) 

mp mp mp 6 mp 

where we have allowed only mod-square combinations of the fields in order 
that the U(l)# and G symmetries remain intact. In a number of recent 
analyses, it has been shown that the scalar spectral index n s can obtain the 



WMAP central value in SUSY models of inflation using a minimal [14] , [15 



or non-minimal [16|, [13[ form of the Kahler potential in Eq. (J4J). The value of 
n s is expected to be an important consideration in determining the correct 
model of inflation; however, these models have predicted exceedingly small 
values of r, and will be handily ruled out should Planck measure a sizable 
tensor amplitude. As we will see, both a non-minimal Kahler potential and a 
substantial inflaton soft mass are needed in order to produce large r-values. 

The F-term SUGRA scalar potential is derived from Eqs. ([3]) and (j4j) 
according to the formula 

V F = e K/m P (K^D Z WD Z *W* - 3m p 2 \W\ 2 ^j , (5) 



where Zj 6 {s, 0, </>,•■■ } are the scalar components of the superfields S, $, 

d 2 K 



$, and where we have defined 



dzidz* 



dW ,dK, 



D Z ;W = (D Z ,W)\ 
Along the D-flat direction, |0| = \<f>\ and Eq. (J5J) takes the form 

V F = n 2 M A (l - J-^ + + ■ ■ ■ V« M 2 (2 (\s\ 2 - M 2 ) + ■ ■ • ) +• 

\ mp 2 mp J v v ' ' 



3 



where 

7 

7s = 1 - 7j K s + 2/t| - 3k SS - (6) 

For simplicity, we will take Ks, k-ss > throughout our calculations. Suitable 
initial conditions ensure that the system of fields evolves along the valley of 
local minima located at |0| = \<f>\ = 0, |s| > s c = M. (For a detailed 
analysis of initial conditions in hybrid inflation, see Ref. [17].) As long as 
the couplings remain perturbative and \s\ < nip (such that SUGRA correc- 
tions remain under control), the potential is dominated by the constant term 
Vq = hz 2 M A . Thus V > and SUSY is broken during inflation, and radia- 
tive corrections and soft SUSY-breaking contributions must be taken into 
account. It has previously been shown in the literature that each of these 
additional contributions can play an important role in the predictions of the 
model Gi- 



lt is convenient to reparametrize the field in terms of a dimensionless 
quantity, x = \s\/M. After including the radiative and soft corrections, and 
retaining terms up to order |s| 4 from the SUGRA contribution, the scalar 
potential during inflation becomes 



V ~ K M 1 - K S x 2 + 75 — + 



X 

\mp J ' J \mp J 2 




where 



. . Iff a (X 4 -1) 2l X 2 + l , K 2 M 2 X 2 \ , s 

F(x) = - 4 ^ + l) In + 2x 2 In + 2 In - 3 J (8) 

encompasses the radiative corrections and 

a = 2 |2 - A\ cos[arg s + arg(2 - A)} (9) 

is the effective coefficient of the linear soft term having a coupling (A — 2) 
in the Lagrangian. The coefficient is the dimensionality of the gauge 
multiplets 0, under G, Q is the renormalization scale, and we take a grav- 
itino mass of 7713/2 ~ 1 TeV. The form of the soft SUSY-breaking terms (i.e. 
the last two terms appearing in Eq. (CO) ) is derived from a gravity-mediated 
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SUSY-breaking scheme, and both of the coefficients a and M| may be either 



positive or negative. Ref. [15| has shown that a negative soft mass-squared 
for the inflaton having a magnitude at intermediate scales can readily lead to 
good agreement with the WMAP central value of n s . For large tensor modes, 
we will see that an intermediate-scale soft mass remains important, yet in the 
present case we will have M| > as in split supersymmetry models 20] . We 



take a at constant values, which can be achieved via an appropriate choice 



of initial conditions [13(. The values we consider will be discrete, with the 
intention of representing the possible sign choices for a. 

The vast majority of models in the literature treat the inflationary dy- 
namics using the slow-roll approximation, in which cosmological quantities 
may be expanded in powers of the slow-roll parameters 

ml fV\ 



AM 2 \ V 
ml fV" 



2M 2 \ V 



(11) 



where primes denote a derivative with respect to x. Within the slow-roll 
approximation (i.e. for e, \rj\, £ 2 <C 1), inflation lasts for a number of e- 
foldings given by 



m \ 2 r Xo ( v 



N ° - 2 L w dx ' (13) 

where x e parametrizes the field value at the end of inflation, and a subscript 
'0' corresponds to the pivot scale ko = 0.002 Mpc -1 crossing the horizon. The 
value of x e is fixed either by the breakdown of the slow roll approximation, 
or by a 'waterfall' destabilization occurring at the value x c — 1 if the slow 
roll approximation holds. To leading order, the scalar spectral index, tensor- 
to-scalar ratio, spectral running, and primordial curvature perturbation may 



5 



(respectively) be written as 



dn 



dink 



n s ~ 1 — 6e + 2?7, (14) 
r ~ 16e, (15) 

i 16er) - 24e 2 - 2£ 2 , (16) 



A 2 - _M^(V1_) (1 7) 

For comparison with WMAP 7-year measurements [2l|, we will evaluate these 
functions at the pivot scale xq. 

In order to construct a viable model it is, of course, necessary to en- 
sure that Eqs. f fT3|) - f|T7|) yield values in good keeping with the most recent 
experimental results. However, if all of the independent parameters in the po- 
tential ([7]) are allowed to vary simultaneously, an analytical solution quickly 
becomes intractable. (This is particularly true if we wish for some of the 
derived quantities to fall within a range of good values, rather than taking 
on the central value alone.) One possibility for circumventing this issue is 
to make various simplifying assumptions to relate or eliminate some of the 
parameters. It is, however, difficult to see how this approach may lead to a 
tensor-to-scalar ratio that is multiple orders of magnitude larger than what 
has been predicted in other treatments of this model. In the next section, 
we will outline the approach we have employed to search for large r-values. 



Seeking Large Tensor Modes 

Our goal is to determine whether there exists a region of the parameter 
space {k, M, ks, k S s, a, Ms, %o} specifying the potential in Eq. (J7J), which can 
lead to large primordial gravity waves. For this, it is necessary to explore 
large regions of this many-dimensional space, and it is advantageous to place 
only the most conservative constraints in order to allow for a thorough inves- 
tigation. While sophisticated computational techniques would likely provide 
a great deal of useful insight, we have determined that a simpler approach 
would suit our criteria quite well. 

In order to explore the parameter space, we have employed a 'brute force' 
random generation of points in this space. For each of the parameters in the 
potential, we have chosen a range within which to generate values randomly. 
For parameters which were expected to vary over multiple orders of magni- 
tude, we have instead chosen to randomly generate the base-10 logarithm of 
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parameter 




type 


quantity 




K 

M/mp 


[10- 4 ,5] 

[io^io- 1 ] 


log 
log 


n s 


[0.920,1.016] 
= 0.968 ± 4a 


Ms/mp 


[io- 8 , 10~ 4 ] 


log 


A' 2 


[2.21,2.65] x 10" 9 




[KT 5 ,3] 


log 




= 2.43 x 10~ 9 ± 2a 




[10- 5 ,3] 


log 


r 


< 1 


a 
x 


{-1,0,1} 


linear 


N 


[50, 60] 



Table 1: Ranges specified for the fundamental parameters in Eq. ([7]), and constraints 
placed on derived quantities. Note that a was considered at discrete values, and xq can take 
on any value between the waterfall point and the Planck scale. Central values and standard 
deviations for measured quantities are in reference to the WMAP 7-year analysis 21] . 



the parameter (e.g. if we allow k- values on the interval [10 -4 , 1], it is per- 
haps more useful to generate logfc on the interval [—4, 0]). Table [1] shows the 
ranges of the fundamental and derived quantities that will correspond to the 
figures in this letter. 

Each set of values {«, M, Kg, uss, a , Ms, Xq} represents a point in param- 
eter space, and specifies a unique potential via Eq. ([7]). After generating 
these values for a given point, functions derived from the potential can be 
calculated with ease. Before considering whether the point fits our desired 
constraints, we check a handful of basic properties. For instance, the physical 
inflaton mass-squared 

+ M 2 S (18) 

is required to be positive at the global minimum of the potential, (s, <fi = 0) — > 
(0,M). (While such considerations are built into the physics of the model, 
there is no guarantee that the points randomly generated in parameter space 
will lead to physical results, and so this type of constraint must be imposed 
by hand.) Also, since we expect the region of interest to coincide with large 
values of the field, the linear term will be suppressed unless a is very large. 
Thus while we retain this term in our calculations, it can effectively be ignored 
for qualitative (and most quantitative) considerations. 

In addition, it is necessary to place constraints in order to ensure that the 
fields evolve to the (SUSY) global minimum. In particular, whenever negative 



m 



inf 



K M 



2-k s ( 
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contributions to the potential are present, there exists a possibility of a false 
(metastable) vacuum developing. If the inflaton is trapped in such a vacuum 
state, inflation may last for a very long time and SUSY will remain broken. 
In order to ensure that the inflaton evolves to the SUSY minimum, we require 
that the potential be essentially monotonic over the interval x e [l,xo]. 

An immediate consequence of the assumption that the potential is mono- 
tonic over a given range is a limitation on the behavior of the polynomial 
terms in V. (The term arising from radiative corrections yields only a sub- 
tle contribution to V, especially at large values of x.) To facilitate this 
discussion, we may rewrite the potential in the form 



where we have defined y = \s\/mp, with y < 1 in order that SUGRA cor- 
rections remain under control. We expect large r-values to be obtained for 
2/o — I? where the quadratic and quartic terms may come into competition 
depending on the size of their coefficients; the quadratic term dominates un- 
less 1 7s | is rather large. Eq. (jSJ) tells us that 75 may in principle be positive 
or negative. For k s , k S s of order unity or smaller, a large positive 75 can 
only be obtained for small values of both couplings; then small values of 
M| are required to obtain an overall quadratic coefficient that is negative. 
For 7s < 0, the quadratic coefficient must be positive, and the combination 
|/?s + 3kss cannot be too large lest the quartic term begin to dominate for 
x < xq. (A similar potential having positive quadratic and negative quartic 
contributions was treated in Ref. |22j.) The case where both the quadratic 
and quartic terms are positive typically leads to n s > 1, owing to r\ > in 



Continuing in this same vein, it will be informative to derive an approxi- 
mate analytical expression for r. As we have already argued, for large 7/0 we 
may suppress the radiative correction and soft linear terms. Using Eq. (fl9|) . 
we may write 



Approximating V(x ) ~ Vq and using Eq. ( ITTj) allows us to eliminate the 




(19) 



Eq. dH]). 
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quadratic coefficient in favor of 77: 




• y [-2751/ 2 + v] ■ 



(20) 



It is worth pausing at this intermediate step to draw a few conclusions. We 
see that, in the large-r regime, the sign of V is determined by a competition of 
two terms. (It turns out that, despite its required smallness during inflation, 
77 is of the same order as the other term here and thus cannot be omitted.) 
Then the requirement that the potential be monotonic up until x = xq leads 
to 



For T] > 0, 75 can be positive, but the largeness of y and necessity of \rj\ <C 1 
demand that [75 1 be very small. On the other hand, 77 < requires 75 < 0, 
and 1 75 1 is not limited by the size of \r)\. That being said, if |7s| is too large 
a local maximum may develop for x significantly removed from the Planck 
scale. Since we require a monotonic potential for x < Xq, this drives Xq to 
smaller values and radiative corrections can no longer be reliably suppressed. 
More stringently, smaller y is expected to lead to smaller values of r. Thus 
we expect large r-values to arise primarily for 75 negative and of modest 
magnitude. 

With these considerations in mind, we may write down an approximate 
formula for r using Eqs. (TTUl [TBI |2"U|) : 



The value of r] at the beginning of inflation is determined in a complicated way 
by the variation of many parameters. It may be useful to seek a substitution 
for a parameter whose value is better known. Indeed, using Eqs. (TH f [T5 | [21]) . 
we find 



It turns out that, after satisfying the full set of constraints, r + is quite large; 
indeed, in most cases this solution violates the slow roll approximation that 
we have assumed in this derivation, as can be seen in Fig. [Q In contrast, r_ 
retains its validity and is an excellent fit in the region where r is large. 



-2 ls y 2 + 77 > 0. 



r ~ IQyl [-2 ls y 2 + 77] 



(21) 



r±-A 1 + 3y 2 (l - n s ) + 12 7 ^ 4 ±yjl + 6^(1 - n a ) + 24 75 |/ 4 . (22) 
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Figure 1: Approximate values r± (see Eq. (J22J)) vs. calculated values r of the tensor-to- 
scalar ratio. The two solutions r+ and r_ are represented by green crosses and purple 
horizontal lines, respectively, with a black dashed line denoting equality with r. We see 
that r_ is the appropriate choice, yielding an excellent approximation for large tensor 
modes. The results displayed here correspond to G = U(l), but the conclusion holds 
equally well in the flipped SU(5) case. 



Results and Discussion 

In order to perform the calculations, it is necessary to make a choice of 
gauge group G, which determines the size of the radiative correction term 
in the potential. We will chiefly discuss the simplest choice of G, namely 
U(l) (i.e. M = 1), which may be identified with a B — L gauge symmetry. 
Some degree of caution must be employed, however, due to the formation 
of cosmic strings upon the breaking of G at the end of inflation. Indeed, 
since Gfi ~ (M/mp) 2 |23( and since large r implies large vacuum energy 
Vq^ 4 = a/kM, we immediately expect high string tension (larger than the 
current bound Gjj, < (2-7) x 10~ 7 2^) to accompany large tensor modes. 



To address this, one may employ a modified framework such as the so-called 
shifted hybrid inflation 0, If, for the purpose of inflating away the strings. 
Although we do not pursue this option, we expect the shifted inflation case 
to yield qualitatively similar results to those presented here. An alternative 
is to make a different choice for G, such as flipped SU(5) j25[, which we will 
discuss at the end of this section. 
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The results of our numerical calculations using G = U(l) are presented 
in Fig. HJ Panel (a) displays the behavior of the tensor-to-scalar ratio with 
respect to the spectral index. We see that these parameters are essentially 
uncorrelated in our model; in particular, large values of r can be obtained for 
essentially any n s value within the 4<r range we have explored. It is especially 
worth pointing out that this model has no difficulty in generating a spectral 
index near the WMAP7 central value n s ~ 0.968. This is due largely to 
negative values of V" at the start of inflation. Referring to Eq. (T141) and 
noting that (typically) e <C \r)\, we see that r\ < will drive the spectral 
index down to red-tilted values (n s < 1). In this regard, these SUSY hybrid 



models resemble models of hilltop inflation [26l . |22 



The first panel of Fig. [5] also serves to define the color- and symbol-coding 
used in the other panels. It is useful, for example, to see that points having n s 
within la of the central value (red X's) cover essentially the same regions as 
the 4<T points (blue crosses) in many slices of the parameter spaced] Perhaps 
even more telling is the location of the points with large r- values, denoted by 
green horizontal lines and for which we have chosen r > 10 -4 as an arbitrary 
cutoff. 

In panel (b) , we show the dependence of the spectral running dn s / d In k on 
r. While the largest r- values can lead to an alarmingly large running, this is 
not necessarily the case; it is also possible to generate large tensor modes for 
much smaller magnitudes of the spectral running. It is not a priori obvious 
which case is more likely. As our aim is to explore the parameter space, we 
simply note that there exists a possibility for observable r with dn s /dlnk in 
keeping with current experimental results. 

Panel (c) shows an overall inverse correlation between k and M, which 
is expected from the estimate V 1,/4 = ^[kM ~ 10 15 ~ 16 GeV. The well-known 
slow-roll formula V(x ) 1/4 ~ 3.3 x 10 16 GeV-r 1 / 4 , with V ~ V(x ), dictates 
that r should obtain its maximal values when one of n, M is large. Based 



upon the reasoning put forth in Refs. [13|, [14J, [15| , one may anticipate large 



tensor modes to occur for large values of k. However, we see here that large 
r is obtained quite readily for moderate values of k, corresponding instead 



1 As can be seen in Table the 4ct range was imposed by taking (central value ±4cr), as 
an initial cut to restrict the results to interesting values. In contrast, what we refer to here 
as the la range corresponds to the region bounded by the la curves in the (n s ,r) plane 



as presented in Ref. [21|. These points were extracted upon analysis of the 4<r results, at 



which stage a reduced number of points made it realistic to use a more precise approach. 
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Results for G = U(l) 





(e) (f) 

Figure 2: Results of the numerical calculations for the U(l) case. As demonstrated in 
panel (a), the points are color- and symbol-coded to represent regions of interest. Blue 
crosses represent all points generated according to Table [T] red X's correspond to a la- 
range of n s ; and green horizontal lines signify points with 'large' r-values (> 10~ 4 ). 



12 



to M-values larger than anticipated, and even approaching the Planck scale. 

From panel (d), it appears that r and Ms are somewhat correlated, es- 
pecially for r > 10~ 3 . This can be understood by noting that the quadratic 
coefficient must be large enough to overcome the influence of the (negative) 
quartic term and stabilize the potential, as discussed in the previous sec- 
tion. In order for this to be true as yo (and thus r) increases, Ms must 
also increase substantially, and becomes quite sizable as the field amplitude 
approaches the Planck scale. In these models, we find Ms/mp < 10 -5 (or 
Ms < 10 13-14 GeV), much larger than the soft mass values utilized in many 
implementations of supersymmetric models. Indeed, these intermediate mass 



scales are reminiscent of split-SUSY models |20j, in which the scalar soft 
masses may lie many orders of magnitude above the fermionic soft masses, 
which remain around TeV-scale. 

Panels (e) and (f ) shed light on the behavior of parameters related to the 
non-minimal Kahler potential. Panel (e) supports our previous claim that 75 
should be negative (and fairly small) at the largest values of r. In panel (f), 
we see that large r- values are mainly obtained for large Kss, in particular 
K ss ^ 1/3 for which 75 < at small k s . It is also interesting to point out 
that if both k,s and kss are sufficiently small, n s cannot lie within the la 
range and r > 10 -4 is not possible. In other words, the use of a non-minimal 
Kahler potential is critical for obtaining large tensor modes, particularly with 
favorable n s . 

As can be seen in Fig. [T] and panels (a), (b), (d), (e) and (g) of Fig. |5J our 
calculations for G = U(l) have yielded values of the tensor-to-scalar ratio 
up to r ~ 0.03, which essentially coincides with the estimated threshold 
for detection by the current Planck satellite experiment. If detected, tensor 
fluctuations in the CMB can give much information about the energy scale 
of inflation as well as gravitational wave amplitudes in the early universe. A 
measurement by Planck lying close to this threshold could serve to support 
the validity of SUSY hybrid inflation models, yet a detection of r significantly 
higher than 0.03 would be difficult to reconcile with the predictions of the 
model that we have presented here. 

The so-called '^-problem' has materialized as a stumbling block to many 
SUSY models of inflation, and as such, it is germane to discuss this issue in 
the context of the current model. In SUGRA hybrid inflation with a minimal 
Kahler potential, this problem is eliminated by a fortuitous cancellation of 
the troublesome mass term with a term arising from the superpotential 
Generalizing to a non-minimal form of the Kahler potential introduces a new 
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mass term for the inflaton, whose contribution cannot be suppressed for k$ 
non-negligible. In our numerical calculations, we have required r\ < 1 for 
long enough to produce 50-60 e-foldings of inflation, i.e. the //-problem is 
not present in the numerical results we have presented in this letter. This 
is accomplished by a delicate cancellation of terms in V"(xq), namely those 
arising from the quadratic and quartic terms in the potential. 

It should also be noted that the behavior described in this letter is, at 
least in part, a consequence of our choice of signs for various parameters. For 
example, our simplification ks, k-ss > leads to the restriction Mj > as 
well as affecting the behavior of the quartic coefficient 75-. Accordingly, it is 
possible that there exist other regions of parameter space leading to sizable 
r which correspond to different sign arrangements of the couplings. We have 
tested this in limited capacity, and while no such regions were revealed, a 
more thorough investigation is needed. It would be particularly interesting 
if such sign changes could result in large tensor modes for Ms ~ 1 TeV. We 
find that this may be possible for the choice k>s < 0, uss > 0, while other 
sign choices come into conflict with our previous assumptions. In this case, 
however, it is unclear whether large r-values can result from perturbative 
values of the couplings. Also, the drastic reduction in the magnitude of Ms 
may lead to the soft linear term playing a more important role and thus 
altering the results. 

The dynamics of the reheating phase can be an important consideration in 
any inflationary model. This is especially true in models that include super- 
gravity, where overproduction of TeV-scale gravitinos can spoil the success 
of Big Bang nucleosynthesis. In our case, large r can readily correspond to 
a large inflaton mass, which in turn can lead to a high gravitino density via 
thermal processes (see Refs. 0, Ifl for reviews) and through non-thermal de- 
cays of the inflaton [27j]. Possible mechanisms for suppressing the gravitino 
abundance include a period of thermal inflation and subsequent dilution by 



entropy production [28]. In this case, we may consider an additional sym- 
metry whose breaking occurs spontaneously at an intermediate scale, while 
the corresponding phase transition is delayed to low energy scales. Such a 
mechanism will also lead to suppression of the baryon asymmetry 28], such 
that the typical leptogenesis schemes [29j (for non-thermal leptogenesis, see 



Ref. [30]) will be insufficient; then, the lepton asymmetry should be over- 
generated in order to emerge at the appropriate abundance after dilution. 
This may be the case, for example, in models of resonant leptogenesis (3l| . 
While the Lyth bound given in Eq. §2§ applies to a generic inflation model, 
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a more stringent version 



may be used in a large class of models |26|, . The additional assumption 
here is that the relative slope |V/V|, and hence also e, is monotonically 
increasing. While this condition is true in many of the inflation models in 
the literature, it need not hold in general; indeed, one may readily verify that 
Eq. ( )23|) cannot yield values r ~ 0.03 for sub-Planckian values of the field. 
Then the model currently being considered must violate this assumption in 
order to obtain the largest (and observationally most interesting) values of 
r that we have generated^] This essentially corresponds to a sign change of 
V" along the inflationary trajectory, causing e to decrease over some region 
of x. (For a class of models exploiting a non-monotonic variation of e, see 
Ref. O-) 

In hybrid models, inflation may end via waterfall (induced by the dynam- 
ics of 0, 0) and it is not necessary for e to increase to unit order. However, 
a drastic decrease in e will cause A^ ~ V/e to increase dramatically, and 
one may begin to worry about the overproduction of primordial black holes 
(PBHs). Conventional wisdom dictates that the amplitude of the curvature 
spectrum A^ should not exceed unity at the end of inflation (when excess 
PBHs can no longer be inflated away). We have checked our numerical results 
and verified that A^ < 10~ 3 , so that PBH overproduction does not concern 
us here. This may be qualitatively understood by considering Eq. ( TT3|) in 
the form Nq ~ j x ° dxj \fe. If e is very small over a significant range of x, we 
see that N will tend to drift outside of the usual 50-60 window unless the 
range of integration ( narrowed. But this is proportional to the A0 

contained in Eqs. ([2]) and (1231) . and a marked decrease in this range will drive 
r to smaller values as well, possibly into a region compatible with Eq. f l23|) 
for which no region of decreasing e is needed. Incidentally, an overdensity 
of PBHs could also be diluted by entropy production in the same way as 
gravitinos, as described above. 

Finally, we turn our attention back to the issue of topological defect pro- 
duction. Since the standard hybrid scenario offers no mechanism to suppress 



2 We thank David Lyth for drawing our attention to this feature of our model, and its 
possible implications including the following discussion on overproduction of primordial 
black holes. 
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Results for G = SU(5) x U(l) 
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(e) (f) 

Figure 3: Results of the numerical calculations for the flipped SU(5) case. The points are 
color- and symbol-coded in the same fashion as Fig. [2j While some parameters experience 
marginal shifts versus the U(l) case, the overall conclusions remain unchanged. 
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the density of defects produced (at the end of inflation) by the breaking of 
G, it may be advantageous to choose a gauge group whose breaking does 
not produce topological defects. To this end, we also consider the so-called 
'flipped SU(5)' (G = SU(5) x U(l) x ) gauge group [25|], corresponding to 
Af = 10. As an added bonus, it has recently been suggested that SUSY 
hybrid inflation models employing a flipped SU(5) GUT may lead to good 
predictions of the proton lifetime [15|. As can be seen in Fig. |3J the change 
M — 1 — > 10 produces small shifts in some parameters, but does not substan- 
tially affect our overall predictions. Indeed, even the quantitative elements 
of the U(l) discussion above remain largely intact. In particular, the flipped 
SU(5) case again leads to r-values up to ~ 0.03, within potential reach of 
the Planck satellite. 

Summary 

Supersymmetric models of hybrid inflation, highly motivated by connec- 
tions with mainstream particle physics, remain in good agreement with ex- 
periment. We have shown that the typical predictions of the tensor-to-scalar 
ratio r <C 1 can be assuaged via the use of a non-minimal Kahler potential 
in conjunction with a sizable soft inflaton mass. We obtain values r < 0.03, 
which may be measurable by the current Planck satellite experiment. This 
conclusion is made possible when the inflaton soft mass Ms is quite large and 
the couplings associated with higher-order terms in the Kahler potential are 
allowed to have a significant magnitude. (Note that there may exist solutions 
with Ms ~ 1 TeV if ks < 0, but this possibility remains to be explored.) 
In addition, large tensor modes are obtained mainly for a concave-downward 
potential at the start of inflation reminiscent of hilltop inflation, which drives 
the spectral index n s toward red-tilted values. The running of the spectral 
index, | dn s /d In k |, can be < 0.01 in agreement with the latest experimental 
results. 
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